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Abstract 

A remarkable multiple analogue of the Stirling numbers of the first and sec- 
ond kind was recently constructed by the author. Certain summation iden- 
tities, and related properties of this family of multiple special numbers are 
investigated in the present paper. 
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1. Introduction 

Stirling numbers of the first and second kind are studied and their prop- 
erties are investigated extensively in number theory and combinatorics. One 
dimensional generalizations of these numbers have also been subject of inter- 
est. An important class of generalizations is their one parameter g-extensions. 
Many have made significant contributions to such g-extensions investigating 
their properties and applications. We will give references to some of these 
important work in Section [3] below. 

In a recent paper, the author took a major step and constructed an el- 
egant multiple gt-generalization of Stirling numbers of the first and second 
kind, besides sequences of other special numbers including multiple bino- 
mial, Fibonacci, Bernoulli, Catalan, and Bell numbers [16j . In this paper, 
we focus on multiple Stirling numbers of both kinds, and give interesting new 
identities satisfied by them. 
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The multiple generalizations developed in [16[ are given in terms of the 
gt-binomial coefficients constructed in the same paper. Its definition may be 
written in the general form as 

W ffl t W V l<T<j<n I J 

where // is a partition of at most n parts, z G C n and g, t G C. The w;^ 
function that enters the definition is a limiting case of the BC n well-poised 
symmetric rational Macdonald function W\. Note that this definition makes 
sense even when \x is not an integer partition, but is a vector \i G C n . 

2. Background 

The basic g-Pochhammer symbol (a; q) a may be defined formally for com- 
plex parameters q, a G C as 

(a) Q = (a; g) Q := / a ^°° (1) 

where the infinite product (a; g)^ is defined by (a; g)^ := n^o^ - a ^)- Note 
that when a = m is a positive integer, the definition reduces to the finite 
product (a; q) m = YYk=o 0- ~ a Q k )- An elliptic analogue is defined by 
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(a;q,p) m := l[0(aq m ) (2) 

fc=0 

where a G C, m is a positive integer, and the normalized elliptic function 
6{x) is given by 

9{x) = 9(x;p) := (x;p) 00 (p/x;p) 00 (3) 

for x,p G C with |p| < 1. The definition is extended to negative m by setting 
(a;q,p) m = l/(ag m ; q,p)_ m . It is clear that when p = 0, the elliptic (a;q,p) m 
reduces to the basic (trigonometric) g-Pochhammer symbol ([T]). 



For any partition A = (Ai, . . . , A n ) and t G C, define [39 



(a) x = (a;q,p,t) x :=Y[(at 1 l ;q,p)x r (4) 
fe=i 
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Note that when A = (Ai) = Ai is a single part partition, then (a;q,p, t) x = 
(a; q,p)x 1 = (o)ai- F° r brevity of notation, we also use 



(a^ . . .,a k ) x = {a u ...,a k ; q,P,t)x ■= (ai)\ ■ ■ ■ (a k ) x - (5) 



Recall that we use V to denote [15j the space of infinite lower-triangular ma- 
trices whose entries are rational functions over the field F = C(q,p,t,r,a,b) 
which are indexed by partitions with respect to the partial inclusion ordering 
C defined by 

H C A ^ < A^ Vi > 1. (6) 

The condition that a matrix u G V is lower triangular with respect to C can 
be stated in the form 

u\f, = 0, when \l £ A. (7) 
The multiplication operation defined by 

(uv)x^ := UxisV Ufl (8) 

fiCuCX 

for matrices u, v G V makes V into an algebra over F. 

2.1. Well-poised Macdonald functions 

The construction of our multiple Stirling numbers involves the elliptic 
well-poised Macdonald functions Wx/^ on BC n 15] . These remarkable fami- 
lies of symmetric rational functions are first introduced in the author's Ph.D. 
thesis 13J in the basic (trigonometric) case, and later in 14| in the more gen- 
eral elliptic form. 

Let A = (Ai, . . . , A n ) and fi = (/ii, . . . , fi n ) be partitions of at most n parts 
for a positive integer n such that the skew partition A/// is a horizontal strip; 



i.e. Ai > /ii > A2 > jU2 > • • • A n > /i n > A n+ i = /U n +i = 0. Following [15 
define 



H x/ ^q,p,t,b):= I 

l<i<j<n 



(g/ , i -^- ltJ --i ) ^_ i _ Aj(g A i +A jt 3-j-i 6 )^_ i _ Aj 
( gW - w _ 1 +l t ,-- i _l) w _ i _ A ^ g A i +A J +l i 2-i-i 6 ) w _ i _^ 
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and 



w ( ,v „ , (a: 1 ,ax) x (qbx/t,qb/{axt)) fl 

W x/ ,(x;q,p,t,a,b) := ff V/l (g,p,f,&) ■ (x - 1>ax) (g6g>g6/M)A 



Hi 0(6*1-*) (M- 2i ),, + A 1+1 ' 1 ) 



i=l 



where g, p, *, x,a,b G C. Note that Wa/^(^; g, p, *, a, &) vanishes unless A/// is 
a horizontal strip. The function Wa/ m (?/, Zi, . . . , ^; g,p, t, a, b) is extended to 
t + 1 variables y,zi,...,2/GC through the following recursion formula 

Wx/^y, zi, z 2 ,..., zt\ q,p, t, a, b) 

= w V»(vt~'\ Q,P, t, at2e > btt ) W v fr(z u ...,z e ; q,p, t, a, b). (11) 

2.2. The Limiting Cases 

The Macdonald functions W\ are essentially equivalent to BC n abelian 
functions constructed independently in 35]. The limiting cases defined above 
are closely related to the Macdonald polynomials 3l| , and interpolation Mac- 
donald polynomials [341 ] . 

The following limiting case of the basic (the p = case of the elliptic) 
W functions will be used in our constructions below. The existence of these 
limits can be seen from (p = case of) the definition ffTUj) . the recursion 
formula (11 ip and the limit rule 

lim a M (x/a)^ = (_l)M x H r »00 g »(>*') (1 2 ) 

where = Ya=i^ and = E"=i(* - and n (^') = E*=i Ca)- 

Denote Hx/n(q, t, b) = Hx/^q, 0, *, b), and for i 6 C define 



W\/Jx;q,t) := lim (s |A| lim Wwja;; g, *, a, as) 

s— >oo V a— >0 

= (.^-W+N^a'h^^^^^Ia (13) 



„-l 

( ;r 

The recurrence formula for u>a/ m function turns out to be 



w x/(t {y, z; q,t) = J2 t m - H) w x/ „(yr e ; q, t) w v/fl (z; q, t) (14) 



Similarly, for x G C define the dual function 
wx/^{x;q,t) := lim ( lim Wx/^x; q, t, a, as) ) 

= r n(A) + H + n M ^ (g;t) MA (15) 

The recurrence formula for the dual wx/^x; q, t) may be written as 

tix/fiiv, z \ q,t) = ^2 ^/^(y*"'; q, t) w\/^(z; q, t) (16) 

for y G C and z G C € . We now recall some old, and derive some new basic 
properties of the w function and its dual, and their connections. 

Corollary 1. Let \x be an n-part partition, and x = (x±, . . . , x n ) G C n . 

(1) The and its dual are flipped versions of one another. That is, 

w^x, q, t) = q-^r^-^-^w^l/x, 1/q, 1/t) 

(2) The limit lim 5 _>i w^xt 6 ^ 1 '; q, t) exists. For the particular case when x = X 
is a partition, we use the notation 

«V(gV<»>; l,t) := lim(l - g)-^,^;?,*) (17) 

q-¥l 

Proof. Both properties follow, by direct calculation, from the definition (j!3p 
of Wx/fi, the recurrence relation ffT4|) for wx/^, and limit formula f lT2|) . and 
the flip formula 

x^(x-\q,t), = (-l^q^'k-^frq-^r 1 ), (18) 

The proof also uses the result that, in the limit 

KnH x/ll (t,q) = @? (19) 
ml 

where k is the maximum of the list k = max{Ai — A2, A2 — A3, . . . , A n _i — A n }, 
and m is the maximum of m — max{Ai — fi±, fit — A 2 , A 2 — 1^2, ■ ■ ■ , Hn-i ~ A n }. 
If the second list has a negative number (which means that A/// is not a 
horizontal strip), then Hx/^,(t,q) = 0. 
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Remark 1. We will need the following properties from 16] in what follows. 

(3) If z = xt s ( n \ for some x G C, we get 

Wll {xt 5 ^;q,t) = (-l^xMr^q-M-^ix- 1 )^ JJ ^ (20) 

l<i<j<n ^ 

which, after flipping g and t and using the flip rule (fTSj) . may be written as 
^(xt^);g,t) = g-^(a;;g- 1 ,r 1 ), J] (21) 

(4) The vanishing property of W functions implies that 

w fl (q x t 5 ;q,t)=0 (22) 
when fj, <2 A, where C denotes the partial inclusion ordering. 

(5) Let A be an n-part partition with A n ^ and < k < X n for some integer 
k, and let ) G C n . It was shown in [16 | that 

n 

w- k (x;q,t) = q- nk H(q 1 - k x l ) k (23) 



8=1 



where k = (k, . . . ,k) & C n . 

2.3. qt-Binomial Coefficients 

The multiple Stirling numbers we develop in this paper are closely con- 
nected with binomial coefficients just as in the one dimensional case. Recall 



that gt-binomial coefficient is defined by [16 



Definition 1. Let z — (xi, . . . , x n ) G C™ and fi be an n-part partition. Then 
the qt-binomial coefficient is defined by 

where q, t G C. It should be noted that this definition makes sense even for 
\i G C n . 
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Setting t = q a and sending q — >■ 1 yields the multiple ordinary a-binomial 
coefficients. For n = 1, the definition reduces to that of the one dimensional 
g-binomial coefficients 



(?), 



(q)n-k(q)k 



(25) 



which are also known as the Gaussian polynomials. These are studied exten- 
sively in the literature including but not limited to the works in 0, 0, [23[ 0, 



25|, 122, |24j, [12 



For the special case z — x = (x, . . . , x) 6 C n , we get 

£2n(/i)+(l-n)M 



q.t 



n 

l<j<j'<n 



(26) 



Note that with this definition we can write the terminating gt-binomial 
theorem 16l in the form 



fj,CX ^ 



x 



(27) 



q.t 



We now recall another important extension [16J that generalize the one 
dimensional g-bracket to multiple gt-bracket as follows. 



Definition 2. Let z = 

[z, sip := [z,s,n,q,t]^ 
1 

g"" " " ' 



l X\ , . . . , x r 



\ Then 



n (i 

i=i < k 



-qt 



n 

l<i<j<n 



Wfl (sqH 5 ^;q,t) (28) 



zs called the partition fi shifted qt-number or the qt-bracket. Note that the 
definition combines a multiplicative variable s, and an exponential variable 
z. Depending on the application we often set z = and write (s)^ = [0, s] M; 
or set s = 1 and write [z]^ = [z, 1] M . In the special case when fi = 1, we get 



[z,l,n,q,t]i 



1 — q Xi t r 



(29) 



which is a multiple analogue of the classical q-bracket. 
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The definition [z]^ may also be written as 



l<i<jr'<n 



which reduces to fl28|) when /i = 1. 

Note also that when 2 = (ar, . . . , x) — x e C n for a single variable x G C, 
the /i-shifted qt- number may be written as 



n{(rz^}(*-i»- i .'- i ).=n{¥^ 1 } (3l) 



This definition reduces to the classical g-bracket in the one variable case. 

Note that (x; 1/q, with the reciprocals of g and t, corresponds to a 

multiple basic gt-analogue of the falling factorial x R := x(x— 1) ■ ■ • (x— (n— 1)) 
as opposed to the rising factorial or the Pochhammer symbol. 

3. Multiple basic and ordinary qri-Stirling numbers 

In this section we review the definition and fundamental properties of the 
multiple Stirling numbers of the first and second kind indexed by partitions 



16l | . The classical Stirling numbers of the first kind are defined to be the 



coefficients of the falling factorial x n in the expansion 



x n = nl 



^ ' k=0 



The g-analogue of these numbers are defined in [5| and their properties are 



studies in [33|, |26|, [29|, |41| and the works of others. 

First, we recall the definition of the multiple gt-Stirling numbers general- 
izing the one dimensional g-analogues. 

Definition 3. For an n-part partition X, the qt-Stirling numbers of the first 
kind S\(X, /j,) are defined by 

mi 

[x] x = Y,Q' niX ' ] t 2nipy{n ' lM si{X,^) J]>i ( 33 ) 

/iCA i=l 
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where x G C, and x % = {x, . . . , x, q, . . . , q] G C n with fi[ copies of x for the 
dual partition fi' . That is, 

iFi=n { L q Jn-l (34) 

i=i i=\ v y 7 

Likewise, the qt-Stirling numbers of the second kind S2(A,/i) are defined by 

1=1 ^ ^ 7 j=l ^tCA 



We now revise an explicit formula for the gt-Stirling numbers given in [16 
generalizing one dimensional g-analogues as follows. 

Theorem 2. For n-part partitions v and fi, an explicit formula for the qt- 
Stirling numbers of first and second kind si(A,/i) and S2(A,/i) are given by 

q n(u') t -2n(p)+2(n-l)\fi\ 

Sl (u, fi) = Sl (u, fi, q, t) = — ^— - f^q, t) 

■ ^ M (z/,A)t (1 - n)|A| w M (g- A r 5 W;l,l/t) (36) 

fiCXCu 



and 



s 2 {v,n) = s 2 {v,ii,q,t) : 



-nO')£2rtO)+(l-«)M 



n?=i(i-9t n " i )" i "' i< 

• J](-l)l A lrW«; A (g^);l,t)/ A (g,t)t;(A, / u,g,t) (37) 

where and are as defined above in CorollaryU\ and f(fi), w(A,/i) and 
v(X, /i) are given by 

' («)«-*♦, TT f S 



/o^-n^sr n ^dp-}* 1 n(". -"«)! (38) 

i=l ^ ^ l<i<j<n < ^ >V*-M > i=l 



l<i<J<n ^ ' W W 
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and 



■ ^((?¥ (n) ;g, £) = (-l)Mg^')t-n(M) / A ) (40) 



Proof. We write the IV-Jackson sum 15 in the form 



W x (x;q,p,t,at- 2n ,bt- n ) 

— 1\ f I -.Lj.—i—i+V, 



(a) aCm" 1 *-"" 1 ) a tt f (^^-a, (g&r 4 - J+1 ^ 



n 



(qbs-H-i) x (qbt«s/a) x { (P^) Xi _ x .(qbt-^) Xi+Xj 

V (bs-H- n ),(qbt n /a), " ( (l-fe-^-V^) 2t _ 2 
4t (^- 1 ) M (as- 1 t-"- 1 ) At All (l-fort 1 ^) ig J 



n 

l<t<j'<n 



■ W^qH^-q, t, bt 1 - 2 * 1 , bs~H- n ) ■ W^xs; q, t, as" 2 *" 2 ", bs'H' 71 ) (41) 

Set b = ar in this identity and send a — > 0, multiply both sides by (rt n )' A l 
and send r — >■ oo using the limit rule (lT2~j) to get 

l<t<j<n ^ ^ Mi i J 
(_1)M g 2M t nQx) g nQxQ J (,, %, 

'2-^ (n+n-l\ 11 



■ ( lim lim ^(g¥ (n) ; q, t, art 1 - 2 * 1 , ars" 1 *'")) ■ w„(xs; g, t)s~ H (42) 



\r— too a— >0 



Simplify this in the one variable case z = xt 6 ^ using ( 12T1) . and send s — )■ oo 
to get 

(x;q-\t- 1 ) x = J2 u ^^)^ 1 ( 4 3) 
jitCA 

with the definition of u(X,/i) above. Similarly, apply the shifts a — > as 2 , 
b —7- bs and x — > x/s in ( ]4T1) at the beginning, and follow the same steps 
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except send s to to get 



EX A ' v) 0;e~\*~% (44) 



where v(X,fi) is defined as in the theorem. It is clear, by a change of basis 
argument, that 

S„\ = u(v,\)v(\,n) (45) 

Note that the left hand side (jHJ) does not depend on q or t. Now, flip the 
parameters q — > l/q,t — > 1/t, take limit q — > 1, and multiply and divide the 
summand by YYi=i V(l — qt n ~ % Y i to get 



mi 



fiCA i=l i=l 

Multiply and divide the summand in ( l4"5j) by £( n-1 )M, substitute ( l4T)j) into 
for (xt™ --1 )'^, multiply both sides of this latter identity by fllLi 1/(1 — <?£ n 
to get 

II { (i _ gt n-^ ) r1 )" = E ( II (1 _ qtn-iy, 

2=1 v v ' y fiCv \i=l v ' 



Ml 



■ ^ «( I /,A)r^ 1 >l A l(Bmt;(A,A*,l/g,l/t)))n^> ( 47 ) 

/jCACi/ / i=l 

Multiplying and dividing the summand now by )f~ 2n (' i )+( n ~ 1 )H gives 

/j,Cv i=l ^ y ' 

where Si(v, /i) is as defined in the theorem. A similar sequence of calculations 
show that 



r (1 - xt n ' l ) Vi _ ^ / " (1 - qt n - 1 )^ 
qt 71 '^ ~ ^ (1 - qt 71 -*)"* 



(limn(z/,A,l/g,l/t))t("- 1 )l A ^(A, / i,g,t)) (x), (49) 
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Multiplying and dividing the summand now by gr»0Of2»(")+(i-»»)M 

gives the 

explicit formula for the gt-Stirling numbers of the second kind 

II ix~-aX = S 5" ( "' , '- 2 " M+( "- I,IH ^(". f ) ■ to, (50) 

i=l ^ ^ ' /xCk 

Finally, substituting x q x completes the proof. We conclude by simplifying 
the flips and limits that entered the formulas above. 
It follows immediately from ( Il8p that, if 

h(fi,q,t):= J] / Jg^pW | and g(», q, t) := {qt n ~ l ), (51) 

l<i<j<n *■ ^ 

then 

M/i,l/g,l/t) = t 2n( ^ {n - 1)M ^(/x,g,t) (52) 

and 

9 (pL,l/qA/t) = (-l)Mq-M^)l>M-WMgfa, q ,t) (53) 
Thus, flipping the parameters give 

f_l)H H+n(M')W^)-(i-i)M 
«(A, ft 1/,, = 

■ n ( fa^-'r" }»,(«v>;«.«) (54) 



and 



v(\,ii,l/q,l/t) 

Multiply and divide both by (1 — g) Atl , and pass the limit as g — )■ 1 to get 
\imu(\,/i,l/q,l/t) 



1 (i n_i L- 11 I (t j_i_1 ) u -u- 

1 v l<i<j<n k V 7Ml ^ 

j^i+i 

n-1 



• Km (i - <?r II rr^ — ^(9^ (B) ; i. ( 56 ) 
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where is defined as above. Using the limit rule ( 1T2j) . direct calculations 
give that 

( 1 n_1 1 \ n_1 

lim (i - qY 1 ^- n 7^ — = - ^ 

Hence 

limw(A,/i,l/g,l/t) = (-l)Wt n ^-^Wt5 M (gVW;l,t)/(u,g,t) (58) 
Similarly, 

InWA,//, 1/g, 1/t) = t^-^w^q-H- 5 ^; 1, 1/t) /(//,<?,*) (59) 
which completes the proof. 

Remark 2. The immediate properties of the Stirling numbers established 
in [l6|] are listed as follows: 

(a) The multiple Stirling numbers Si(u, fi) and s 2 (^, fi) admit explicit com- 
binatorial formulas which are derived in the Theorem above. 

(b) These explicit formulas reduce to those for the g-Stirling numbers given 
by Kim in [29J for n = 1. Moreover, sending q — > 1 in that case yields clas- 
sical Stirling numbers of both types. 

(c) The matrix m with entries = si(A, fi) is invertible in the sense of V 
algebra, and its inverse is given by mTi = s 2 (A,/i). More precisely, we have 

5 uX = si(v, A) s 2 (A,yu) = ^ s 2 (z/, A) si(A,/i) (60) 

which follows immediately from the inversion relation (I45p . 

(d) Similar to the one dimensional case for the g-Stirling numbers, we have 

Sl (A,A) = s 2 (A,A) = l (61) 
for an arbitrary n-part partition A. 

(e) Setting t = q a and sending q — > 1 gives the multiple ordinary a-Stirling 
numbers of the first and second kind. 
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4. Additional multiple Stirling number identities 

We derive some additional new properties of the multiple Stirling num- 
bers in this section. 

(1) First note that lim x ^ (x)^ = YYi=i^-/(^ ~~ qt 11 ' 1 )^ follows readily from 
the formula fl3Tj) written in the gt-angle brackets (x)^ as 

(S) = ff J L_l ( x; q-^r 1 ), = 17 ( ^^M (62 ) 

In the multiple case, setting x = in (jJSJ) and (1501 respectively gives 

n n 
i=l v H ' \iC.v i=l v y ' 

and 

TT- V —— = y q n ^h- 2n ^ n - 1 ^s 2 (u,f,)T\- (64) 

These appear to be new identities, even in the one dimensional case n = 1. 

(2) Note that, for an n-part partition v with v n ^ 0, the bracket {x) u has 
roots at x = t x ~iq m i for j = 1, . . . , n, and rrij = 0, . . . , i/j- — 1. The limit 
bracket J^ftt =1 (l — xt n ~*) M4 /(l — gt™ -1 )^ has roots at x = t 1- -' for j = 1, . . . , n. 
Therefore, if we set x = t x ~iq m i (for some rrij < Vj) in (1501) we get 

n 1 1 n ' m 'i+ 1 — j+n— i\u% 

= £ JJ (lzj_J ) ,,(„, „) (65 , 

/x£f i=l ^ y 7 

where the summation is over all partitions // C i>, that is all partitions /jCi/ 
such that /ij < m^. This inequality follows from the vanishing property of 
the w functions (|2"2"|) . 

In the particular case, setting x = q in (T4"gj) gives 

= ^ g -"(^)^W-(r»-i)W Sl ( A) ^ ( 66 ) 
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which is an analogue of Ylk=o s i( n > A;) = in the classical case. 

In another special case, setting x = (i.e., rrij = 0) in and (1511 
would amount to vanishing of all brackets (x)^ except the ones in whose index 
the j-th part (thus all parts j + 1, . . . , n after j) are 0. That is, the brackets 
will be nonzero only for partitions such as /i = . . . , fij-i, 0, . . . , 0), and 
others will vanish. This is particularly interesting, for the substitution x = 
t 1_J the limit brackets niLi(l ~~ xt 11 ^ 1 )^ / (1 — qt 71 ^ 1 )^ also vanish except for 
partitions /i = (/ii, . . . , [i n -j, 0, . . . , 0). 

In general, setting x = t x ~i in (HBI) and (!50|) respectively gives 

where the sum is over all partitions fi C v such that fi = (fix, . . . , fij-i, 0, . . . , 0) 
Likewise, 

= ? "^') r 2nW+(n-l)|,| S2 ^ ) ^ ( f l-^ (6g) 

where the sum is over all partitions C v such that /i = (jjli, . . . , 0, . . . , 0) 
The particular cases when j = 1 in (1671) and j = n in (|68|) show that the 
multiple Stirling numbers vanish when jj, — = (0, . . . , 0) G C n as in the 
classical case. That is, 

Sl (A,0) = s 2 (A,0) = 
for any n-part partition A with A n ^ 0. 

(3) Recall that, Si(n, m) = — s 2 (n, m) = — (™) when n—m = 1 for the classical 
Stirling numbers. Similarly, if the index partitions satisfy |A| — |A| = 1, we 
have that 

Si(A,A) = -s 2 (A, A) 

exactly as in the one dimensional case. 

The proof follows easily from the inversion fl60l) relation, and the observa- 
tion that there are only two partitions between A and A under the inclusion 
ordering, namely the two partitions themselves. That is, 

= 5 AA= Sl(X,(l)8 2 (jl,X) (69) 

AC^CA 

which implies that si(A, A) s 2 (A, A) = — si(A, A) s 2 (A, A). That the diagonal 
entries of both type of multiple gt-Stirling numbers are 1 by (l6Tj) is now 
enough to conclude. 
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5. Conclusion 

We have derived several interesting summation identities for the multi- 
ple gt-Stirling numbers in the present paper. We will write down additional 
properties such as the recurrence relations they satisfy, their explicit evalua- 
tions in certain special cases, their combinatorial interpretation, generating 
functions, and connections to other families of multiple special numbers in 
an upcoming article. The Stirling numbers have interesting connections to 
various branches in mathematics such as the one expressed in the classical 
Dobinski's formula. Such relations will also be formulated in that paper. 
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